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A scheme for retrieving quantum information stored in collective atomic spin systems onto optical 
pulses is presented. Two off-resonant light pulses cross the atomic medium in two orthogonal 
directions and are interferometrically recombined in such a way that one of the outputs carries most 
of the information stored in the medium. In contrast to previous schemes our approach requires 
neither multiple passes through the medium nor feedback on the light after passing the sample 
which makes the scheme very efficient. The price for that is some added noise which is however 
small enough for the method to beat the classical limits. 
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I. INTRODUCTION 

The development of complex quantum communication 
networks requires the interface between light and atomic 
carriers of quantum information. The photons are ideal 
for transmitting the quantum information over long dis- 
tances. The atoms, on the other hand, are optimally 
suited for local storage of quantum states in the internal 
ground-state coherences which can exhibit very long life- 
time. Ideally, the atoms should form a quantum memory 
for light where the quantum state of the optical beam 
can be stored, possibly processed, and retrieved later on 
in a controlled manner. The quantum memory is an es- 
sential ingredient of quantum repeaters 0, 0] that allow 
for long-distance quantum communication over realistic 
noisy quantum channels and it is also necessary for scal- 
able all-optical quantum computing as proposed by Knill, 
Laflamme and Milburn 

Several methods for exchanging the quantum states of 
light and atoms have been considered in the literature 
such as storing the state of a single photon in a single 
atom j4| or coupling the light to an ensemble of atoms 
HUHlllllaliallll- The latter approach seems 
to be particularly advantageous because it is not neces- 
sary to employ a cavity and work in the strong coupling 
regime. Instead sufficiently strong interaction is achieved 
due to the collective enhancement effect [Tmi. An al- 
ternative to direct light-matter state exchange is to first 
prepare an entangled state of atoms and light and then 
teleport the quantum state of light onto atoms. The en- 
tanglement between a photon and a single trapped ion 
has been reported T^l and also entanglement between a 
photon and an ensemble of atoms has been observed |0 . 
In related experiments, correlated photon pairs emitted 
by an atomic ensemble were detected 0, 0, . 

The storage of a quantum state of light mode in the 
ensemble of atoms has been recently demonstrated ex- 
perimentally by exploiting the off-resonant interac- 



tion of light with a cloud of cesium atoms held in a 
glass cell at room temperature 12, 18, 19^, 29, 21] . In 
that experiment, weak coherent states were stored for a 
few milliseconds and the fidelity of the stored states ex- 
ceeded the best fidelity achievable by a classical measure- 
and-prepare procedure. The storage protocol consisted 
of sending the light beam through the atomic ensemble, 
performing a suitable measurement on the outgoing light 
and applying an appropriate feedback on atoms with the 
help of a magnetic field. This procedure, however, cannot 
be directly inverted to retrieve the quantum state from 
atoms back onto light due to the current technical lim- 
itations. The light pulses used in the experiment need 
to be about 1 ms long (i.e. 300 km) and the memory 
readout would require storing somewhere the first pulse 
which interacted with the atomic ensemble while a sec- 
ond auxiliary light beam is used to measure the atoms. 
The required long duration of the pulses also prevents one 
to employ measurement-free retrieval methods based on 
multiple passages of the light pulse through the atomic 
sample (22i.i23i.i24j|. 

A possible solution to these problems has been sug- 
gested by using two- and four-pass protocols with the 
light beam crossing the medium simultaneously from 
two perpendicular directions 25J . Although high fidelity 
readout can theoretically be achieved in these protocols 
they each suffer from deficiencies which can decrease their 
general applicability. The two-pass protocol is asymmet- 
ric with respect to the two conjugate quadratures and a 
unity gain readout is difficult to achieve since it would re- 
quire a complicated temporal profile of the coupling light 
beam. This makes the two-pass protocol unsuitable for 
readout of general input states. With the four-pass pro- 
tocol perfect readout can be achieved in principle, how- 
ever, since it involves so many passages even a small light 
refiection loss at each crossing can seriously decrease the 
achievable fidelities. This point will be investigated quan- 
titatively for our protocol below. 
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In this paper we suggest an alternative protocol for 
read-out of the atomic quantum memory which requires 
only a single passage of the light through the atoms and 
is therefore rather robust with respect to losses on the 
glass cell walls. The simultaneous transfer of both atomic 
quadratures xa and pA onto light is achieved by sending 
two light beams perpendicularly to each other through 
the memory unit, see Fig. ^ The scheme can be ex- 
plained on the simplified model with a single atomic en- 
semble as in Fig. however, to avoid technical noises, 
experiments are usually done using two cells in external 
magnetic field (see Fig. ^p), observing the phenomena on 
the sidebands of the hght pulses. As it turns out, these 
two approaches yield similar results, however, the num- 
ber of crosses of the beams through the cell walls plays 
a crucial role in increasing the noise. 



The rest of the paper is organized as follows. In Section 
II we derive the formula for the output light quadratures 
and the readout fidelity. In Section III we take into ac- 
count the effect of losses and show how they influence the 
read-out fidelity in both the single-cell and double-cell 
schemes. We will also demonstrate that the readout can 
be improved by squeezing the input light beams entering 
the cell. Finally, Section IV contains brief conclusions. 




FIG. 1: Scheme for the single-passage read-out of the atomic 
quantum memory. Two light beams L and M simultaneously 
pass through one (a) or two (b) atomic ensembles. The wave 
plate (WP) swaps the quadratures x and p in one of the 
beams. Afterwords, the beams are (re)combined on a bal- 
anced beam splitter BS. 



II. READOUT PROTOCOL 
A. Single-cell scheme 

Let us first deal with the single-cell scheme with no 
magnetic field. The beam L propagates along the z axis 
and the Hamiltonian describing its coupling with the 
atomic ensemble is = aS^Jz and the second beam 
M propagates along the y axis and interacts with the en- 



semble according to Hm 
Hamiltonian is 



aS^ Jy so that the complete 



H = a {S^Jz + Si' J, 



(1) 



Here S denotes the Stokes vector characterizing the 
polarization properties of a light beam, J stands for 
the collective atomic spin operator and a is the cou- 
pling constant ITM l20j |. The orientation of the 
Stokes vectors is chosen such that S3 describes the he- 
licity and Si describes linear polarization, with posi- 
tive 5*1 corresponding to a;-polarized light. The com- 
ponents of S and J satisfy the angular momentum- 
type commutation relations, [Jj,Jk\ — iejkiJi with 
j,k,le {x,y,z} and [Sj^{t), S^it')] = ie,kiS{t - t')St {t), 
and [Sfit),Sf{t')] = ieju5{t ~ t')S^'{t) with 3,k,l e 
{1,2,3}. 

Both beams M and L contain strong x polarized com- 
ponent such that the components of the Stokes vectors 
Si{t) and Si\t) can be approximated by c-numbers. We 
assume that the Si components are constant throughout 
the duration of the pulses and both beams have the same 
intensity, 

5f (t) = S^'it) ^^^Si, (2) 

where is the total number of photons in the strong 
component of each beam and T is the duration of the 
pulse. Similarly, the atoms are polarized along the x axis 
and the x component of the collective atomic spin attains 
a macroscopic value and the operator can be replaced 
with a c-number, « (Jx)- 

In the Heisenberg picture, the atomic quadratures 
evolve according to 



^ = ~a{Jx)Sr{t)- (3) 



The relations between input and output y and z com- 
ponents of the Stokes vectors of beams L and M read 

El, 



aSiJy{t), 



(4) 

The Heisenberg equations (O can be easily solved and 
yield 

Jy{t) ^ Jy{Q) + a{Jx) f Si{T)dT, 
Jo 
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Jo 



(5) 



We define the quadratures of tlie light modes L and M 
in the usual way as integrals over the whole pulse, 



/ 2 f 

XL.M = \ Hr S2''^Ht)dt, 

V J-^L Jo 



(6) 



We a lso d efine the ato mic q uadrature operators xa = 
Jyl \J {Jx)i PA — Jzl \J {Jx)- All these operators sat- 
isfy the canonical commutation relations [xj,pfc] = iSjk- 
Upon combining the equations lO and © and 

changing the order of integration we obtain the expres- 



sion for x™', 



dt\l- 



qM 
^3 



it) 



(7) 



where we introduced the effective coupling constant k — 
a\J Nl{Jx) /2. We can see that the description in terms 
of the light modes ^ is incomplete because a mode with 
non-constant temporal profile 1 — t/T couples to x™*- 
The two mode profiles we are dealing with are 



hit) 



/2(i)=V3 1 



(8) 



with the normalization \f'j{t)\dt = T. The two mode 
functions ^ are not orthogonal since fi{t)f2it)dt — 
Ta/3/2. We can decompose /2(i) as follows, 



/2W = ^/lW 



1 



fsit), 



(9) 
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(11) 



Note that the information about the xa and pA 
quadratures is stored in two different modes L and M. 
In order to obtain a single mode that would contain an 
approximate version of the state stored originally in the 
atomic memory, we have to combine the modes L and 
M. To accomplish this, we suggest placing a wave plate 
to one of the beams, say M, to switch the quadratures 
xm Pm, Pm —xm, and to let these modes inter- 
fere on a balanced beam splitter, as shown in Fig. ^ In 
one of the output arms we finally obtain a mode with 
quadratures 
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71^ 



1 



p ■ 



V2 



ip' 



,out 



^m')- (12) 



On inserting the formulas Hll|) into (|12|l . we find 
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^(Y-I)Pl 



(13) 



If we choose k = v^, then the quadratures are transfered 
from the atomic memory to the light mode with unit gain. 



1 



1 



Pa 



V2 ^ V6 



Pm, 



P 



-x^-—x,,~—p^. 



(14) 



where fzit) = V3(l — 2t/T) and the functions /i(t) and 
f^it) are now orthogonal We define quadrature operators 
corresponding to the modal profile /3(i), 



xl,m 



Pl,m 



L Jo 



2t 



^ 1-=^ S. 



T 
2t 



-iL.M 



-iL,M 



{t)dt, 



{t)dt. (10) 



This is the main result of this paper. 

The memory read-out is not perfect due to the noise 
of the light quadrature operators in Eq. (|14|l . To be 
specific, let us consider a readout of a coherent state la) 
stored in the atomic ensemble. If the light modes are 
initially in vacuum states, then the added noise is the 
same for both quadratures and can be expressed in terms 
of the mean photon number n of thermal light with the 
same quadrature fluctuations 



One can easily check that the quadratures (|l()|l satisfy 
canonical commutation relations. 

Using the above definitions we can express the output 
quadrature operator xj"* as a linear combination of the 
input light and atomic quadrature operators. This pro- 
cedure can be repeated for other quadratures and after 
some algebra we finally get: 



KP2 



-Pm 



2V3 



P^M, 



pT 



PL, 



The fidelity of the memory read-out is in this case 

1 



1 + n 



(15) 



(16) 



and for n = 1/3 we get T = 3/4. This value significantly 
exceeds the best fidelity J^mcas = 1/2 achievable by mea- 
sure and prepare procedures, which in '2^ was shown to 
be the optimal classical storage or retrieval procedure. 
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Note that the readout of the memory is imperfect be- 
cause a part of the information on the state remains in 
the memory. In particular, for k = V2 the final atomic 
quadratures read 



L I 



„out 

Pa 



Pa 



V2p] 



M- 



(17) 



The distribution of the quantum information between the 
quantum memory and the output light beam suggests 
that the readout protocol is related to the (asymmet- 
ric) quantum cloning. Indeed, in the limit of infinitely 
squeezed quadratures pL = Pm = 0, the formulas ifT^ 
and l(T7|) describe exactly the two clones produced by 
the optimal Gaussian asymmetric 1 — > 2 cloning ma- 
chine for coherent states The anti-clone (or, more 
precisely, the complex conjugated clone) is carried by 
the light beam in the auxiliary output port of the beam 
splitter, whose quadrature operators can be expressed as 



Pa 



1 
1 



1 



(18) 



atomic sample 2 is then described by 

= aiS^jf^ + Sl'jf^) + ^Jf\ (20) 

The information is encoded in the sum- or difference- 
quadratures of the collective atomic spins 



XA± 
PA± 



1 ( 



1 Ml) -U (2 
-= [p\'±p\ 



V2 



(21) 
(22) 



The input-output relations can be derived in a similar 
manner as for the singic-ccU scheme by solving the linear 
Heisenberg equations of motion. In the limit f^T ^ 1 
which is satisfied in the experiments (fJT w 300) one 
thus gets 
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B. Two-cell scheme 

To avoid technical noises, one can work with a two-cell 
scheme as shown in Fig. 9]. In this model we assume 
magnetic field B perpendicular to the two beams such 
that the fields in the two cells have opposite orientation. 
The atomic spins then rotate around the a:-axis with the 
Larmor frequency fl (X B, the direction of rotation being 
opposite in the two samples. The atoms interact with 
light on the sidebands with frequencies loq ± f2, where loq 
is the carrier frequency. The quadratures of cosine and 
sine modes at sideband frequency ft ^ l/T are defined 
as follows. 



■'^L,M 



rP 

Pl,m 



^L,M 



Pl,m 



2 
2 
2 



S^-^'{t) cos{iU)dt, 



5*3 '^"^(t) cos{nt)dt, 



S^'^'^it) sm{nt)dt, 



S^-^^t) sm{nt)dt. 



The quadratures xf^'fj and p]^'f^ corresponding to modes 

with temporal profile V3{1 — 2t/T) can be defined in a 
similar way. The interaction between light and atomic 
sample 1 can be described by means of the Hamiltonian 



(19) 



where the term proportional to arises due to the ap- 
plied magnetic field. The interaction between light and 
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for the light beams, and for the atomic samples 
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(23) 



(24) 



So as to have the atomic quadratures in one beam, we 
place a quarter-wave plate in beam M, thus switching 



the quadratures x'^^j^ — > p^f and p^^^ 



-x.f , and 



"-M I'm I'm ^ -^M 

combine the L and M beams on a beam splitter. In one 

of the outputs we then find the new quadratures 
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(25) 



in J_ C,in 
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^M 



C,in 

Pl 



" =P^^'", (26) 



2V6 
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~S,in 
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V2 
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^SAn 
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~S,in 

Pl ■ 



(27) 



(28) 



If we choose k — \/2, the input-output relations for the 
hght are 







1 C.in 

x/2'^^ 


1 ~C,in 
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1 ~C,in 






— ~^A- 


, 1 s.in 
H l=Xr 

V2 ^ 


1 -S,in 




p' 


= -p'X^ 


1 S.in 


1 -s.in 


(29) 



which is completely analogous to the single-cell results in 
Eq. p4|) . In this case the atomic sum- and difference- 
quadratures are written into the cosine and sine side- 
bands of one of the output beams. 



III. LOSSES AT THE CELL WALLS 
A. Vacuum input modes 

So far we have assumed perfect light transmission 
through the cell walls. In reality, passing through the 
boundaries between different media causes losses which 
result in decreasing the read-out fidelity. To estimate the 
effect of losses, each boundary is modeled by a beam- 
splitter which transforms the input light quadratures as 
= Vl - Ax™ + VAx^""", where A is the effective loss 
coefficient, x™'^ is the vacuum-noise field quadrature and 
analogous transformation applies for the p-quadrature. 

In this subsection we assume the input modes L and 
M in vacuum states. We shall consider read-out with 
unity gain so that the mean values of the quadratures 
are faithfully transferred from the atoms onto light and 
the readout fidelity of coherent states \a) is constant and 



does not depend on a. The field attenuation \/l — A can 
be compensated by increasing the coupling constant k, 
but the extra vacuum noise decreases the quality of the 
read-out. It is also possible to counteract the losses by 
amplifying the output light beam in a phase insensitive 
amplifier with gain G = 2/[k^{1 — A)]. However, this also 
adds some thermal noise. 

After suitable modification of Eqs. Hll|l and we 
find that for A < 2/3 the amplifier does not help and the 
best result is achieved with the choice k = ■\/2/(l — A) 
which leads to the read-out transformation with added 
noise corresponding to mean thermal photon number n 



squeezed input light 
single-cell scheme 




FIG. 2: Dependence of the fidelity of coherent state read- 
out T on the loss coefficient A at the cell surface is shown 
for vacuum input light field and both the single-cell scheme 
(solid line) and the double-cell scheme (dashed line). The dot- 
dashed line depicts the fidelity achieved when using optimally 
squeezed input light in the single-cell scheme. 



and coherent state read-out fidelity J- given by 
1 _ 3(1-^) 



3(1 -A)' 



4- 3A 



(30) 



Note that for small losses one gets the fidelity approxi- 
mately 



3 

4 



3A 



(31) 



On the other hand, for A > 2/3 we find that the opti- 
mal coupling constant is = 2-^/3/(1 — A), the output 
signal needs to be amplified, and the resulting noise and 
fidelity read, 



^3(1-^) 



- 1, 



T=-^'3il-A). (32) 



The dependence of on A can be seen in Fig. [21 (solid 
line). As can be seen, the scheme is rather robust with 
respect to losses so that fidelity above 1/2 is achieved 
with the loss coefficient up to A = 2/3. 

When using the two-cell scheme (Fig. H^), the influ- 
ence of losses is more complicated. To make the discus- 
sion simple, let us consider the case when the mean values 
of the quadratures are retrieved with unit gain without 
amplification of the output beam and we also require that 
the cosine mode of the light beam does not couple to 
the difference-quadrature atomic mode with quadratures 
XA~ and PA-- These conditions uniquely determine the 
values of the coupling constants which must be different 
in the two cells, in particular, we have 



(1-A)3 



K2 



(1-A)- 



(33) 
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The calculation of the read-out fidelity can be most eas- 
ily carried out with the help of the formalism of Gaus- 
sian completely positive maps p^ . Efll 13(1 which al- 
lows us to describe the effects of losses in a particu- 
larly simple way. Assuming the atomic modes as well 
as the light beams are initially in coherent states, the 
joint state of the atoms-|-light system remains Gaussian 
throughout the evolution and can be fully characterized 
by the first and second moments of the quadrature oper- 
ators. Moreover, since some quadratures remain uncou- 
pled during the readout procedure it suffices to consider 
only seven quadratures that can be conveniently collected 
into a vector r = ,pf The 

covariance matrix which captures the quadrature fluc- 
tuations and correlations is defined as follows, jjk = 
{Tjfk + fkfj) — 2{rj){rk). The initial covariance ma- 
trix is the identity matrix, = 1. After the passage 
through the first atomic ensemble, 7 transforms accord- 
ing to 71 — SASijinSf S]l + G, where 



Si 



describes the coupling of light with the atoms and the 
matrices Sa = Vl — ^Is® I2 and G — ^l5©02 account 
for the losses on the cell wall. Here Ij denotes an identity 
matrix with j rows and columns and Oj stands for j x j 
square matrix with all elements equal to zero. After the 
passage through the second atomic sample the covariance 
matrix reads, 

72 = SaS2{Sa7iSI + G)S^Sl + G. 

This formula accounts for the losses on the two walls of 
the second glass cell as well as the interaction with the 
second atomic ensemble, 
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The final interference of the beams L and M on a bal- 
anced beam splitter can be again described by a symplec- 
tic matrix but for our purposes it suffices to evaluate the 
variance of the resulting quadrature operator x'~' which 
is related to the thermal noise added to retrieved state 



via n = {{Ax^)'^) - 1/2. It holds that 
n = 7(72,11 + 72,22 + 272,12) 



After tedious but straightforward algebra one finds the 
added thermal noise and read-out fidelity to be 



1-A- 



3(1-^)3 ' 



3(1-^)3 



4-10^+^^2-3^3 



Thus, for small losses one finds 



4 8 ' 



.(34) 



(35) 



i.e. , the influence of losses is twice as large as in the single- 
cell scheme. Note that the required ratio between the 
two coupling constants, k,2/ki = (1 — A) is naturally 
achieved in practice because k oc V-^l and the intensity 
of the light beam is attenuated by a factor of {1 — AY 
when crossing the two cell walls that separate the first 
and second atomic ensembles (c.f. Fig. 1). The fidelity 
as a function of A is plotted in Fig. |21 (dashed line). As 
can be seen, with large losses the fidelity decreases much 
faster in the two-cell scheme than in the single-cell one. 



B. Squeezed auxiliary modes 

So far we have assumed that the unused input light 
modes were in vacuum states. In principle, the fidelity of 
the transfer can be improved by using light with squeezed 
quadratures xl, Pl, xm, and pM- Let us first con- 
sider the probably experimentally simplest case where we 
squeeze the x quadratures of all modes in both beams L 
and M, irrespective of their temporal profile. Assuming 
that pure squeezed vacuum with the minimum uncer- 
tainty is used, we have 



((A<)2) = ((A4'^)2) 



-2r 



The fidelity reads 



12 



(36) 



(37) 



(38) 



The optimal squeezing that maximizes the fidelity is 
e2'- = 73 and 



6 



6 + V3 



0.776, 



(39) 



which is only very small improvement over the read-out 
with light beams in vacuum states. 

Much better results can be expected when squeezing is 
very selectively distributed among the modes such that 
XL, xm, Pl, and pm are squeezed. In such a case the 
fidelity can approach unity with increasing squeezing, 
provided losses are negligible. Presence of losses at the 
boundaries decreases fidelity (see Fig. \^ . For the sake of 
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simplicity let us consider the single-cell scheme. Optimal 
read-out is achieved when the quadratures pM and pL 
are very strongly squeezed, while the optimal squeezing 
of the modes L and M is finite for < ^ < 1 and the vari- 
ances of the four quadratures can be expressed in terms of 
a single variance F as follows, {{Axl)^) ~ {{Axm)"^) = V 
and ((Api)2) = ((Apm)^) = 1/(4^). We find that for 
A < ^th = 5(3 — ^7/3) it is optimal to choose 



while for A > ^th we get 
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A 



2{1~A)' 



3(2- A) 
1-A 
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In this latter case, it is necessary to amplify the signal 
in the output light mode in a phase insensitive amplifier 
with a gain G — 2/[k^(1 — A)]. The resulting optimal 
fidelity can be obtained from Eq. H16|l where the mean 
number of noise photons reads 



n — < 



A, A>At^. 



2 A, 



I 2-A 
3ir~A) 



(40) 



One finds that for A 1 the fidelity decreases with the 
loss coefficient as w 1 — (4/3)^. 



be very efficient since only a single pass of the beam 
through the medium is needed. Thus, no storage of the 
long (300 km) pulse is necessary and the losses at the 
media boundary are minimized. The price for this oper- 
ational simplicity is some impcrfcctncss of the read-out if 
vacuum input fields are used: providing no losses occur, 
the read-out fidelity of coherent states can approach the 
limiting value 3 /4 which is still substantially higher than 
the measure-and-prepare limit of 1/2. Unit fidelity can 
be approached in our scheme if strongly squeezed fields 
in the relevant input modes are used. We have studied 
the influence of losses at the cell walls, showing that the 
fidelity decreases with losses twice as fast when a two- 
cell scheme is used in comparison to a single-cell scheme. 
This also confirms our motivation to search for simple, 
compact schemes for quantum information exchange be- 
tween field and material carriers. Our results can lead 
to building quantum memory devices suitable for storing 
and straightforward retrieving of information encoded in 
continuous variables. 
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